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Abstract 

Chiral perturbation theory supplemented by the Omnes function 
is employed to study the strength of the isoscalar central nuclear in- 
teraction, Gs, in the chiral limit vs the physical case. A very large 
modification is seen , i.e. ?y s = Gs chirai/Gg physical = 1.37 ± 0.10. 
This large effect is seen to arise dominantly at low energy from the ex- 
tra contributions made by massless pions at energies near the physical 
threshold where the physical spectral function must vanish kinemati- 
cally. The slope away from the chiral limit, ds, is also calculated and 
is correspondingly large. I also explain why this large variation is to 
be expected. 
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1 Introduction 



For the most part we have obtained only a phenomenological description 
of the nuclear force. Whether using meson exchange potentials or effective 
field theory, the parameters of the internucleon interaction are hard to relate 
directly to QCD. However I will show that we know enough to understand 
the quark mass dependence of the nuclear central interaction with reasonable 
control. I will use this to describe the strength of the nuclear interaction 
in the chiral limit. We will see that there exists a strong variation of the 
strength, which nevertheless comes from readily understandable physics. 

The basic framework will use a dispersive representation for the scalar- 
isoscalar interaction. The overall strength will be governed by a contact 
interaction 



The strength of this interaction will be related to a dispersion relation^ |2] 



Here the spectral function in the integrand psit- 1 ) is to be calculated from 
physical intermediate states of two-pion exchange, using chiral amplitudes 
at low energy plus the Omnes function for pion rescattering. We can con- 
trol the mass dependence of these ingredients to a great extent using chiral 
perturbation theory. 

Let us at this stage show the basic result, to be developed more fully 
below. The spectral function for the physical case was developed in Ref. [3] 
where it was shown to reproduce the shape and magnitude expected from 
past experience with scalar exchange. This result, with a small modification 
to be described below, is shown in Fig. 1. Also shown is the result that will 
be developed for the spectral function in the chiral limit. The integral under 
theses curves gives the strength of the scalar interaction 1 . The chiral limit is 
seen to have a greater strength largely because the threshold extends down to 
zero energy and the spectral function develops quickly. As will be discussed 
more fully below, there is some modeling involved in producing these spectral 
functions, and there are some gaps in our understanding of mass effects. 
However, the bulk of the change in the chiral limit comes from relatively 
low energies, where we have better control over the calculation. These low 
energy effects are easy to defend as predictions of chiral perturbation theory. 

lr The negative "peak" in this figure is phenomenologically equivalent to the exchange 
of a "er" resonance 



H contact = G s NNNN + ... 
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Figure 1: The spectral integrand, ps(n)/fi, whose integral determines the 
strength of the scalar-isoscalar interaction. The curve that extends down 
to jjL = (solid line) is that for the chiral limit - the other (dashed line) is 
calculated with the physical pion mass. 



These contact interactions enter modern descriptions of nuclear binding. 
They have been developed in a systematic fashion in the treatment of light 
nuclei. Below I will also discuss the modification in the binding energy in 
heavy nuclei using the contact interactions. We will find a large increase in 
the binding energy per nucleon. 

In chiral perturbation theory, the coupling constants will have an expan- 
sion in the pion mass. The leading terms will be 

G s = Gsv{l + dsml)+F s ml + ... (3) 

Equivalently one sees also the parameter Ds = dsGs defined. I will use my 
results to generate a value for the chiral slope, ds- 

The results indicate rather large shifts as one goes to the chiral limit. In 
particular, I will find 

Gslchiral -. — , n 1 n 
VS, ch = = L37 ± °- 10 

S\phy steal 

0.31 ±0.08 

ds = 2 4 

When unexpectedly large effects are found, it is important to carefully un- 
derstand their origin and assess whether the effects are reasonable. I spend 
considerable effort on this task in this paper and we can in fact understand 
why a large effect is found. 
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In section 2, I describe some properties of effective field theory that make 
it reasonable that such a large mass variation exists. In Section 3, the basic 
framework is displayed. Section 4 discusses the leading chiral amplitudes for 
the spectral function. Section 5 describes the Omnes function and the the 
framework that I use for approximating the phase shifts. Section 6 puts these 
ingredients together and shows the comparison of the physical case to the 
chiral limit. I also discuss the uncertainties of the calculation here. Section 
7 uses this scalar coupling to describe nuclear binding in the chiral limit. 
The slope away from the chiral limit ds is extracted in Section 8. Finally 
in section 9 I briefly summarize the results and discuss the relation of this 
work to previous treatments of the nuclear force in the chiral limit [3J IE]. 

2 Effective field theory and low energy effects 

There is a dichotomy between the way that we describe the philosophy of 
effective field theory and and the way that we implement it in practice. The 
difference is important for the understanding of the present problem. Let 
me introduce these notions through a simpler problem which is by now well 
understood, and which also has an anomalously large mass dependence - 
namely the pion and kaon electromagnetic mass splitting. I then show how 
this philosophy is relevant to the problem of the scalar interaction in nuclei. 

An effective field theory is the description of the low energy limit of 
a theory, with the correct degrees of freedom and interactions for the low 
energy world. As originally explained by WilsonJJj, we describe the low 
energy interaction below some scale A by the propagation of the light degrees 
of freedom. The high energy effects above the scale A are incorporated into 
a set of local Lagrangians. They must be local at energies E « A due to 
the uncertainty principle. We then have a clear way of separating low energy 
effects, which we can readily calculate, from high energy effects, which we 
often cannot calculate but which can readily parameterized. 

However, this is not how we generally proceed in practice. The use of a 
cutoff A to separate low energy from high energy is awkward. Energy cutoffs 
require care in order to not violate Lorentz invariance, gauge invariance 
and / or chiral symmetry. In practice we generally calculate using dimensional 
regularization. However, this scheme does not have a separation between 
low energy and high energy - all scales are integrated over in loop diagrams. 
This is not a problem in principle. Extra or missing contributions can be 
adjusted by a shift in the values of the contact interactions. However, and 
this will be our topic here, unexpectedly large mass effects can occur if large 
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low energy contributions are put into the contact interaction. 
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Figure 2: The photon loop diagram which contributes to the pion and kaon 
electromagnetic mass differences. 

Let us explore this first in the context of the pion electromagnetic mass 
difference. This is described by the chiral Lagrangian 



in usual notation. One consequence of this is the equality of the electromag- 
netic mass splitting of pions and kaons 



which is known as Dashen's theorem [Hj. Normally corrections to this would 
be expected to be of order 25% from the naive dimensional analysis estimate 
of SU(3) breaking due to quark mass differences. However direct calculation 
of the mass differences in chiral-based models yields a 100% violation of 
Dashen's theorem^- The reason why the models are right and the naive di- 
mensional analysis estimate is wrong has to do with the distinction described 
above. 

In the effective field theory the degrees of freedom are pions, kaons and 
photons. In a Wilsonian effective field theory, the quantum effects of the 
light particles in the diagram of Figure 2 should be calculated up to the 
scale A, and only the effects beyond that scale parameterized by the chiral 
Lagrangian with coefficient gEM{^)- Thus in a Wilsonian scheme there 
would be two different contributions: a dynamical contribution of the actual 
7r, K, 7 loop diagrams up to the scale A and a parameter 5em(A) describing 
the physics beyond the scale A. There will also be higher order terms in 
the Lagrangian with extra factors of the pion and kaon masses. While we 
expect only a modest variation of the chiral Lagrangian with quark mass - 



Cem = 9em Tr(QUQU ] ) 



(5) 




(6) 
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this is the real content of Dashen's theorem - there is no such guarantee for 
dynamical effects of the light particles. 

In dimensional regularization the structure is different. The photon loop 
in dimensional regularization has the form (for the pion self energy) 



e 2 m 2 
16?r 2 d - 4 



(7) 



The factor of m 2 is forced by dimensional analysis - there is no other dimen- 
sional factor in the calculation. Because of this factor, the loop does not 
renormalize gEM, but goes into the renormalization of a higher dimension 
term in the lagrangian that has extra factors of m 2 and m 2 K - let us call this 
parameter Qem m 2 - So in dimensional regularization there is no residual 
dynamical loop contribution, and the analysis only involves the parameters 
of the chiral Lagrangian treated at tree level. 

-Q 2 U V -a{Q 2 ) 

.008 




,4 0.6 0.8 1 1.2 1.4 

Q 2 (GeV 2 ) 



Figure 3: The figure shows the vacuum polarization function Q 2 IIv—a(Q 2 ) 
whose integral determines the chiral parameter qem- Notice that the great- 
est contribution comes from low energy. 

It is easy to demonstrate that the greatest contributions to the chiral 
coefficient qem comes from low energy physics, not from high energy. There 
is a dispersive sum rule ^Uj for gEM, or equivalently the pion electromagnetic 
mass difference in the chiral limit, for which data exists. In the chiral limit, 
we have 

g EM = -r- / dQ 2 Q 2 U V _ A (Q 2 ) (8) 
Jo 

where Ily-^ is the difference of the vector and axial-vector vacuum polar- 
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ization functions 2 . In turn there is a spectral sum rule for ILy-A, and the 
input to this sum rule can be obtained from r decay data. The integrand, 
derived from ALEPH data, is shown if Fig. |31 We see directly that most of 
the input comes from low energy. 

Now consider what happens if we adopt the Wilsonian scheme and the 
dynamical low energy contribution is relatively large. This dynamical effect 
can easily (and does) have over a 100% variation when comparing the pion 
with the kaon, because m 2 K ~ 13m 2 and we are in an energy region com- 
parable to the masses. In this case, even if the mass dependence of short 
distance physics parameterized by the chiral Lagrangian is of normal size, 
the kaon mass shift can be very much different from that of the pion. 

The Wilsonian approach is undoubtedly correct in saying that there is 
large SU(3) breaking at low energy. It is readily calculable using just the low 
energy theory. However, when doing chiral perturbation theory in the usual 
fashion, regularized dimensionally, we would miss this effect. Dimensional 
regularization can accommodate this situation, but it is not the outcome 
that we would naively expect. It can be done only at the expense of having 
an anomalously large parameter in the chiral Lagrangian. The mass shift is 
described by Qem an d Qem, m 2 > which are a priori unknown. The latter must 
be taken to be much larger than naive expectation in order to reproduce the 
large difference. 

What has happened is that, in this scheme, important low energy physics 
has been encoded in the parameter of the chiral Lagrangian Qem m 2 - We 
normally expect that the chiral parameters are manifestations of only high 
scale physics but in a scheme like dimensional regularization, that does not 
have a separation between high and low scales, there can also be low energy 
physics encoded in the parameters. 

The lesson from this is that if low energy dynamical effects make a sig- 
nificant contribution to a process there can be a large mass effect. This can 
lead to a larger-than-expected variation in chiral parameters when treated 
in the usual fashion. I will now argue that an analogous effect will occur in 
the parameter Gs describing the strength of the scalar interaction. That is, 
we will see that low energy effects are embedded in this parameter and that 
these effects can have a very large variation. 

Consider now the two pion exchange contributions to the internucleon 
potential. The diagrams and formulas will be given later. In a Wilsonian 
approach we would calculate explicitly the pionic effects up to some separa- 

2 There is an extensive literature on Tlv-A - the figure is drawn from my own work on 
the subiectfTT], 
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tion scale A, and then use a contact interaction Gs(A) to describe the high 
energy effects beyond this scale. We might realistically expect that the mass 
effects in Gs(A) are modest. However, there can be no such expectations 
for the explicit two pion calculation if there are significant contributions 
from energies around the pion mass. These contributions can vary by 100% 
as the pion mass is set to zero. For example, in Figure 1 we see that the 
physical spectral function can only start at \x = 2m n while the chiral limit 
case extends down to zero energy. The two cases do not get close to each 
other until above 400 MeV, and the integrated values of the contribution 
below this energy are significantly different. 

However this Wilsonian approach is usually not carried out in practice. 
The interaction in the scalar channel is parameterized by a contact inter- 
actions (or a sigma potential) and the low energy parts of the two pion 
exchange are not separately evaluated. One calculates nuclear binding to- 
tally with the contact interaction or the potential. The consequence of this 
is that we should expect that the large mass variation of the dynamical 
two-pion exchange has to be included in the contact interaction, which then 
appears as if there were an anomalously large mass dependence of this coef- 
ficient. We see in Figure 1 that the high energy portion has a modest mass 
dependence while the low energy component has a 100% variations between 
the physical case and the chiral limit. If we were following the Wilsonian 
approach, the parameter Gs(A) would have a small mass variation if A were 
chosen above 500 MeV. However, in a usual calculation the mass variation 
of Gs, encoded in the parameter ds, must be taken to be very large in order 
to account for the low energy effect. It is this which allows us to understand 
the large mass effect described in the following sections. 

It should be said that there is now an attempt at what I call the Wilso- 
nian approach. Epelbaum, Glockle and Meifiner jl21 113j have been calculat- 
ing the pionic effects to the nuclear potential in chiral perturbation theory 
using a cutoff in the spectral integral. When supplemented with a cutoff de- 
pendent contact interactions this is the procedure described above. It would 
be good to see this approach extended to the binding of heavy nuclei. It 
would be even preferable, without loss of generality, to use the Omnes for- 
malism developed in the present paper for the low energy amplitude, even if 
a cutoff is used. The Omnes function is a required feature of the full answer 
and tames the growing polynomial behavior of the chiral amplitudes such 
that there is not an excessive contribution from the energy region around 
the cutoff. This procedure is equally general because any mistake that is 
introduced from the region around the cutoff can nevertheless be corrected 
for by an adjustment of the contact term. 



7 



3 Formalism 



Our procedure begins with the dispersion relation derived by Cottingham, 
Vinh Mau and others F° r textbook reviews, see |14l I15j . The scatter- 
ing amplitude for two nucleons obeys an unsubtracted t-channel dispersion 
relation. 

,w s 1 f°° , 2 ImM(s,n 2 ) 

M(s,t) = - dp 2 1 (9) 

hm\ P -t-ie 

The imaginary part of this amplitude is connected to the crossed channel 
NN — > NN with the important intermediate state being that of two pions. 
The overall amplitude is decomposed into partial waves described by their 
spin and isospin quantum numbers. The greatest interest in this paper will 
be on the scalar-isoscalar (J=0, 1=0) channel. By taking the nonrelativistic 
limit and ignoring the energy dependence in the S channel, one can define 
a momentum space potential that depends only on the momentum transfer 
q 2 . Using the scalar isoscalar central potential as an example, let us define 
the corresponding spectral function 

Psip) = ImV s {q = ip) (10) 

In terms of this imaginary part the potential is defined by the dispersion 
relation 

When converted coordinate space this defines an inter-nucleon potential 
depending on the separation r, with a spectral representation 

Vsir) = -^- f°° dfifi e~» r ps(fi) (12) 

For orientation, note that this description would produce the conventional 
sigma exchange potential with the substitution 

Ps(p) = -KsfcSQj? - ml) (13) 

This would recover the classic Yukawa potential in momentum and coordi- 
nate space: 

v ° = 2~i ■ ( u ) 

and 

V s (r) = -A- (15) 



S 



However, the ^-function is not a good representation of the physics in the 
scalar channel. Instead we will use chiral perturbation theory extended to 
higher energy with the Omnes representation to describe the ingredients to 
the scalar channel. 

For all channels except that of one-pion-exchange, the potential is short 
range, and in effective field theory can be represented by a delta function - 
a contact interaction. For example, the scalar and vector isoscalar contact 
interactions would be represented as 

H con tact = GsNNNN + GyN^NN^N + ... (16) 

Since the Fourier transform of a constant is a delta function, the contact 
interactions in coordinate space correspond to a constant in momentum 
space. These contact interactions are therefore given by the momentum 
space potential at q 2 = 0, i.e. for the scalar or vector channel 

2 f°° tin 

G S y = V s ,v(q 2 = 0) = - / PS y{n 2 ) (17) 

ft J(2m 7r ,3m 7T ) f 1 

Higher powers of momenta can be accommodated by derivative contact in- 
teractions. For the exchange of a narrow resonance, the contact interaction 
has the form G% = ±g 2 /m 2 , and I will use this relation for the vector meson 
effect 3 However, for the scalar interaction it is important to have a more 
complete evaluation, as described below. 

4 Chiral amplitudes 

The application of chiral perturbation theory to the nuclear force has an 
extensive literature - for reviews and references, see |131 1171 [T51 119j . The 
pathway using the dispersion relation and the spectral function has been 
pioneered by Kaiser, Meifiner and collaborators 201 1211 1221 1T2] . We will find 
this a very effective description of the nature of the energy expansion for 
the nuclear interaction. 

The low energy behavior of the two pion exchange diagrams have been 
calculated in perturbation theory. The ttNN vertex is proportional to the 
axial charge gA and the two pion vertices are parameterized by low energy 
constants c%, C2, C3, C4 in the chiral Lagrangian. The diagrams are shown 

3 For resonance saturation estimates of these and other nuclear contact interactions, see 

El 
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Figure 4: Two pion exchange diagrams which arise in chiral perturbation 
theory. 



in Fig. 3. The diagrams of Fig 3a, b lead to the spectral function [27 
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The result of [20] was calculated in the heavy baryon limit and I have mod- 
ified their result to include the effects of a finite nucleon mass see 
also|24|. The difference is minor numerically. Diagram 3b is one power 
higher in the energy expansion and involves two factors of the low energy 
constants 01,2,3 The imaginary part of this diagram is |21j 



p|(m) 



32vrFi 



4m 2 

IF 



9{fji — 2m n ) 



(19) 



4 Cl ml + ^(/i 2 - 4m 2 ) + c 3 (/x 2 





2mi 



r 2 



4m 



2\2 



4 Watch out for a sign difference - their momentum space potential has the opposite 
sign of the conventions used in the present paper. 
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Figures 3d,e, the box and cross box diagrams, have an ambiguous status 
in the calculation of the contact interaction. In an effective field theory 
treatment, one includes the ttNN vertex in the low energy Lagrangian and 
treats pion exchange dynamically, including loop processes. This procedure 
would generate these two diagrams explicitly. Therefore in principle, the 
correct treatment is to treat the box and cross box separately from the 
contact interaction. In practice is is not always clear how much of these 
diagrams are calculated dynamically and how much is parameterized. This 
probably accounts for the relatively wide range of values of Gs that one 
can find in the literature. In a relativistic treatment such as applied in 
Walecka-style models which I use in Sec. 8, one would expect that these 
terms should not be included in the contact term. This is because the pion 
and nucleon propagators are treated relativistically, in which case the box 
and crossed box diagrams would be calculated directly and should not be 
double counted by also appearing in the contact term. For this reason we will 
not include them in our calculation of the contact interaction. We are also 
fortunate that, on a purely practical level, these terms are small compared 
to the effects that we do study, at least for the scalar central potential. 
This is demonstrated in 20 and and can also be seen in Fig 3.15 of Ref. 
|15j . As an explicit example, if take the q 2 = limit of the potential from 
these diagrams quoted in [201, one finds that they generate a shift in the 
non-analytic mass correction of size 

AGs = 15 f^L = 0-42 GeV- 2 (20) 

* 1024ttF4M v ' 

where the total result is Gs ~ —400 GeV -2 . For these reasons we will not 
use Fig 3d,e in our calculation of the contact interaction. 



5 Unitarity and the Omnes function 

A key ingredient in the present method is the use of the Omnes function, 
which incorporates the physical effects of pion rescattering^j. This ingre- 
dient is not optional - it is required by unitarity. Moreover the general form 
of the result is well know within the elastic region. The amplitude must 
have the form of a polynomial in the energy times the Omnes function 

nM = ^/*-fiiL] (21) 

7T J S S — /X z 

Here 5 is the tttt scattering phase shift, in our case for the 1=0, J=0 channel. 
Chiral perturbation theory is consistent with this order by order in the 
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energy expansion. Following Ref. it is known how to match this general 
description to the results of chiral perturbation theory by appropriately 
identifying the polynomial. The two pion vertices that appear in Fig. 3a,b,c 
are modified by this Omnes function in a way that is strictly analogous to 
the two pion vertex described in[2fi]. 

Experience has shown that using the lowest order amplitudes supple- 
mented by the Omnes function gives reasonably good results up to 700-800 
MeV. I have elsewhere |3 applied this formalism to the nuclear central po- 
tential using the physical phase shifts matched to the lowest order results 
in order to generate the spectral function. The result is very encouraging 
because it generates the main features known phenomenologically for this 
interaction. It was shown that this give good results for the shape and mag- 
nitude of the spectral function in the scalar channel. In traditional treat- 
ments this shape is described by a broad low mass "<r" particle, although 
it is unlikely that such a resonance is present in the spectrum of QCD. In 
particular, in the treatment of [3| the spectral function is given by 

Ps ^) = pfRen^) + p c s \^)\ 2 (22) 

The effect of the Omnes function is to reshape the spectral function so that 
it mimics a a even though such a resonance is not present in reality. I 
will continue to use this formalism in the present paper. While the use of 
the Omnes function is rigorous and can be matched to chiral perturbation 
theory to any given order, the use of only the lowest order amplitudes in 
this expression is clearly not rigorous except at low energies - it is here that 
the modeling of the spectral function enters. 

We have recorded the chiral amplitudes in the previous section. The 
task now is to understand how these amplitudes vary between the physical 
case and the chiral limit. There will be several ingredients that need to be 
explored. Some mass dependence is explicit in the chiral amplitudes. Other 
dependence is contained in the parameters F w g&. Finally there is mass 
dependence in the Omnes function, or equivalently in the tttt phase shifts. 
This latter feature requires the most work, and I will address it first. 

At low momentum, we have explicit expressions for the tttt scattering am- 
plitudes through the work of Gasser and Leutwvler|27|. At higher energies, 
our best knowledge comes from a treatment that combines chiral symme- 
try, dispersion relations, crossing and experimental data|2Sj by Colangelo, 
Gasser and Leutwyler (CGL). In jSj I used the CGL phase shifts to construct 
the Omnes function. However, because of the reliance on experimental data, 
we don't directly have the ability to vary the quark mass in this analysis. 
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In order to calculate the Omnes function for the chiral limit, we need to in- 
troduce a method to model the high energy behavior of the nir phase shifts 
in such a way that we can vary the pion mass. 

There exists a good and successful method for extending the descrip- 
tion of the chiral amplitudes, often referred to as Pade approximation or 
the inverse amplitude method(IAM) [55]. This provides an analytic approx- 
imation to the scattering amplitudes in a form that is fully satisfies unitary 
and which can be matched order by order to the results of chiral perturba- 
tion theory. Matched to the results at order E , one has the result for the 
1 = 0, J = amplitude is 

t 2 

Too = t 2 + t 4 + ... - — (23) 

12 — t4 

The second term t 4 contains the effects of loop diagrams and hence also has 
imaginary parts. The beauty of the I AM representation is that this simple 
rearrangement allows the amplitude to satisfy unitarity exactly. While the 
method necessarily differs from the full answer beyond the order to which 
has been matched 5 , studies have shown that this representation provides a 
good description of the scattering amplitudes up to reasonably large energies. 
For example, the scalar-isoscalar IAM amplitude is compared to the physical 
result of CGL in Fig [3 The description is quite good up to 700 MeV, after 

5 

80 
60 
40 
20 



0.2 0.4 0.6 0.8 

\i (GeV) 

Figure 5: The irir phase shift from CGL (solid line) and that found from the 
inverse amplitude method (dashed line). The agreement is excellent up to 
about 700 MeV. 

5 For example some two- loop logarithms are missing when the IAM is matched at one 
loop order |3U] 
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Figure 6: The absolute square (a) and the real part (b) of the Omnes func- 
tions derived from the CGL phase shifts (solid line) and those found using 
the inverse amplitude method (dashed line). 



which it falls short. The resulting Omnes function is also compared to that 
derived from the CGL phase shifts in Fig 5. 6 . Again the results are similar, 
except at high energies where the Omnes function is small. 

I will use the IAM phase shifts as an analytic approximation of the true 
amplitudes, keeping in mind the shortcomings at higher energies. The low 
energy behavior of these phase shifts will then be fully rigorous and the high 
energy portion will be approximate. 

Let us immediately look at the phase shift and the Omnes function in 
the chiral limit. The phase is shown in figure [7| in comparison with the IAM 
approximation to the physical case. The salient feature is the threshold 
behavior. The phase shift in the chiral limit clearly extends down to s = 0, 
and since this is the S-wave, the strength turns on relatively quickly. The 
physical case needs to vanish at the physical threshold. These requirements 
naturally yields a larger phase shift for the chiral case throughout much 
of the physical region. At high energies, we see a much smaller difference 
between the chiral limit and the physical case. This is what should be 
expected, as the pion mass should make less of a difference at the higher 
energies. The Omnes function follows directly from the phase shifts and also 
has recognizable features. This is shown in Fig EI 

6 In producing the Omnes function, I had to extend the phase shifts above the fi = 850 
MeV endpoint of the CGL analysis in order that the principle value part of the Omnes 
function integral be well behaved near the upper end. Likewise in the IAM formalism, 
there is a corresponding extension. As long as this extension is smooth it has little effect 
on this calculation. 
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Figure 7: The scalar-isoscalar tttt phase shifts using the IAM in the physical 
case (dashed line) and in the chiral limit (solid line). 

These results are encouraging for the reliability of the method. The place 
where we have the least theoretical control is the upper energy end of the 
energy region. The IAM does not produce a much variation at these energies, 
so that most of the variation that we find then comes from physics at lower 
energies. Moreover, the phase shifts only enter our calculation through the 
Omnes function, and the Omnes function is small at high energy. The true 
Omnes function from the CGL phase shifts is yet smaller, so that the IAM 
does slightly overemphasize the higher energies. However, since there is not 
much mass variation at these energies, this is not a serious flaw. Overall, 
despite the approximate nature of the IAM method, there is no reason to 
doubt the general features of the resulting phase shift and hence the Omnes 
function. 

6 Dispersive chiral analysis 

In this section, the evaluation of the contact interaction will be presented. 
The key features of the scalar channel (i.e. the "sigma") emerge, as reported 
in Moreover, the chiral limit will induce a shift in the scalar coupling that 
again emerges mostly from low energy two-pion exchange. In performing this 
comparison, I have not varied the strange quark mass, keeping it fixed at its 
physical value. 

In addition to the explicit dependence on the pion mass in the chiral 
amplitudes and the Omnes function, there is also implicit dependence on 
the pion mass contained in the parameters qa and F v . The dependence has 
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Figure 8: The absolute square (a) and the real part (b) of the Omnes func- 
tions in the physical case (dashed line) and in the chiral limit (solid line). 
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with 1 4 = 4.4 ± 0.2 and diQ = —1.0 ± 0.7. While the mass dependence of 
the pion decay constant is fairly well constrained, the corresponding result 
in qa is less well understood. This will be a significant component of our 
final error estimate. 

The chiral parameters c\p, 3 are determined from a chiral analysis of pion- 
nucleon scattering [38 , \ 6^. When used without matching to the Omnes func- 
tion at second order in m 2 , as is the case both in the ttN analyses and in this 
paper, there is an ambiguity as to whether the fitted parameter should be 
identified with C3 itself or 03^(2771,7,-) which is the threshold value of the cou- 
pling. Since £l(2m n ) ~ 1.25, this makes a difference in the result. The spec- 
tral function is dominantly determined by the parameter C3, which is however 
not very well determined, with C3 = — 4.7l}'g GeV -2 being the quoted value. 
If we used instead the constraint C3Cl(2m n ) = — 4.7j^; GeV -2 , we would 
choose C3 = — 3.7^Qg. The overall magnitude of Gs does depend sensitively 
on the choice of C3, see 0. However the ratio of Gs in the chiral limit to 
that of the physical limit is almost completely insensitive to C3, as will be 
commented on below. Therefore, I will only quote the ratio as my final 
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result. Although I have explored a wide range of C3 values, the figures were 
produced using C3 = —4.0 GeV -2 , a result intermediate between the two 
choices mentioned above and consistent with both within errors. The other 
parameter choices used were c\ = —0.64 GeV -2 (i.e. c\VL{2m 1T ) = —0.8), 
and C2 = 3.3 GeV -2 . I have explored the sensitivity of the results to these 
parameters, and both the magnitude and the ratio are insensitive to reason- 
able variations. 

At this stage we can put the pieces together. The chiral amplitudes of 
Eq. E3contain both explicit and implicit mass dependence, displayed above. 
The Omnes function has further dependence. The resulting integrand in the 
spectral sum rule for Gs has been previously shown in Fig. 1, both for the 
physical case and for the chiral limit. 

As described in [3~], the shape of the spectral function in the physical case 
provides a reasonable representation of the sigma that traditionally appears 
in the nuclear potential. Moreover, for reasonable values of the chiral pa- 
rameters, in particular C3 which has the greatest impact, the physical value 
of Gs can be reproduced. In this way we can understand the nature of 
"sigma exchange" in the nuclear force. The sigma is not easily understood 
in terms of the quark and gluon degrees of freedom of QCD. It is also not 
a conventional resonance, i.e. a pole in the scattering amplitude close to 
the real axis with the scattering phase passing through 90° near the pole. 
In 7T7T scattering, the phase shift is not large near the sigma and careful 
analysis [33] reveals a pole on the second sheet which is very far from the 
real axis. Other approaches have also described an effective sigma exchange 
without invoking a specific resonance [3""1 1331 131)1 13*7] . In the present case we 
have a very simple set of ingredients, each if which is required within its re- 
spective domain. Even though the amplitude is not conventionally resonant 
at 500-600 MeV, the broad peaking of the spectral function reproduces an 
effect very similar to a sigma resonance. 

Since we have analytic control over the ingredients of this calculation 
we may use this formalism to take the pion mass to zero. The result in 
the chiral limit has a similar high energy behavior and the magnitude there 
is only slightly increased. What is most striking is the enhancement at 
low energy. Much of this is purely kinematic. The threshold in the chiral 
limit extends down to zero energy and the chiral predictions for the spectral 
function develops quickly. 

To calculate the ratio of the scalar strength in the chiral limit to the 
physical limit, one compares the integral under the two curves of Fig. 1. 
The extra weight under the spectral integrand causes the value of Gs to 
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increase significantly. Overall I find 

^Slchiral -. - /oc\ 

Vs, ch = 7^ = 1-37 (25) 

\physical 

This is our primary result. It will be applied to the nuclear interaction in a 
subsequent section. 

Let me attempt to assess the uncertainties in this result. It is clear that 
this calculation is not straight chiral perturbation theory, and hence does 
not share the rigor of that method. I have modeled the moderate energy 
behavior in the dispersive sum rule, and neglected all truly high energy 
effects beyond 1 GeV. We clearly have no control over the effects at very 
high energy. However, there is no indication either phenomenologically or 
theoretically that they are very important in the determination of Gs, nor 
that they would have significant dependence on the pion mass. We can look 
at the low and moderate energy uncertainties in more detail. 

Most of the the shift in the scalar coupling comes from the low energy 
region. The threshold behavior is an unambiguous feature - the chiral limit 
amplitudes must extend down to zero energy and the chiral expansion be- 
comes exact there. Even at moderate energies, the IAM representation of 
the pion phase shifts should be quite good. The greatest uncertainty in this 
region is the lack of understanding of the chiral behavior in the pion cou- 
pling to nucleons, as the leading spectral representation at low energy has 
a factor of g\ . This effect is quantifiable. The uncertainty in the chiral 
parameters d\§ and I4 corresponds to an uncertainty of ±0.07 in our final 
result. Perhaps this uncertainty may be reduced in the future by lattice 
calculations. There could be additional m^/(lGeV) 2 ~ 0.02 corrections to 
the 01,2,3 vertices in addition to that modeled by the Omnes factor, so there 
are clearly other uncertainties at the order of several percent. 

The overall strength of Gs depends most heavily on the chiral parameter 
C3. However, the ratio of the chiral limit to the physical case has only 
a remarkably tiny variation with this parameter. Throughout the whole 
range of C3 considered above, we find the ratio changes by less than 1%. 
Thus the uncertainty due to the magnitude of this parameter is insignificant 
compared to other uncertainties. 

One might worry that the treatment of two-pion exchange as a potential, 
or the use of a contact interaction instead of a potential, would break down 
in the chiral limit because of the massless pions. This could be the case if 
the very long range tail of the potential was significantly modified. Within 
this calculation, there is no evidence of a problem. The spectral function 
also can be used to predict a spatial potential V(r). The physical result and 
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the chiral limit result are shown in Fig. |§J The shape of the chiral result is 
quite similar to the physical limit. 
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Figure 9: The attractive scalar-isoscalar potential in the physical case 
(dashed line) and in the chiral limit (solid line). 

The moderate energy effects beyond the realm of straight chiral pertur- 
bation theory are less quantifiable, but we have reason to believe that they 
are not large. In contrast to the low energy effects, there is no reason to 
expect that these effects are enhanced over the usual expectation of being of 
order m^/A 2 . Indeed we find a reasonably small shift in the spectral func- 
tion at high energy, so that our explicit calculation is consistent with this. 
However, I find it hard to defend that high energy shift as a solid prediction 
of the method. If I were to neglect all shifts above 500 MeV, the main result 
for rjs would be about 5% smaller. 

We can perform a further test by adding some expected higher energy 
effects. One uncertainly that we are able to explore numerically involves 
the energy dependence of the pion vertices. The two pion vertex would 
presumably have some formfactor associated with the leading coupling, C3, 
which would generate a dependence of the vertex on the momentum transfer. 
I have studied this by creating an energy dependent vertex through the 
substitution 



As mentioned in [3J, this modifies the shape of the spectral function some- 
what, and changes the integral that determines Gs- However, the change 
is modest enough that the desired value can still be obtained by chang- 
ing the fit value of C3 within the allowed range. In this case there is a 




(26) 
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modest increase chiral/physical ratio of Gs- For example with n = 1 and 
m = 1200 MeV there is a 29% decrease in the predicted value of Gs at fixed 
C3. However, for the ratio r]s, ch this amounts to a 7% increase leading to 
the result rjs, ch = 1-46. The decrease in Gs and increase in rjs is readily un- 
derstood. The formfactor will decrease the contribution at higher values of 
H, but leave unchanged the threshold effects at small fi. The former feature 
will decrease the physical value of Gs, but the new ingredients enhancing 
the chiral limit will be unchanged. 

Additionally, in the section describing the chiral slope below, I will men- 
tion a puzzle concerning the lack of understanding of the nonanalytic behav- 
ior. While it is numerically small, this is also a feature which must somehow 
be generated from the higher energy portion of the spectral function. 

I would summarize this discussion of the uncertainties by giving a final 
result of 

Vs, ch = ^ s ^ chiral = 1.37 ± 0.10 (27) 

\physical 

7 The vector channel 

Also contributing to the central force is the dispersive channel involving 
three pions in an / = state -i.e. the u exchange channel. In this section, 
I will estimate the modification of this strength in the chiral limit. The 
change is very slight, and the net effect on nuclear binding is so small that 
it is well below even the uncertainty of the result in the scalar channel. The 
reason for this is clear - the 3ir channel is very small at threshold and hence 
the threshold sensitivity that the scalar channel exhibited is absent. The 
chiral modification has the natural size of order m^/(lGeV) 2 . 
In the dispersion relation 

G V = - r - PvV) (28) 

7T J 3 mn V 

the spectral function pvifJ 2 ) is small at low energy. In chiral perturba- 
tion theory it arises first as the cut in a two loop diagram, and hence is 
higher order in the energy expansion. There is also a kinematic suppres- 
sion since the 3 pions are in a total spin-one state. Both theoretically and 
phenomenologically there seems to be little strength in this channel, see e.g. 
|38| . until one reaches the energy of the w(783), which is a narrow resonance 
coupled essentially entirely to three pions. After the cj(783), there is not 
another three pion resonance until well over 1 GeV. It is then a reasonable 
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approximation to take the vector spectral function dominated by the narrow 
resonance a; (783). 

We need to understand how the properties of the w(783) - i.e. its mass 
and coupling - are modified in the chiral limit. We can estimate its mass 
in the chiral limit from the SU(3) breaking in the other vector mesons. 
Replacing one of the light quarks by a strange quark yields the i^T* (890) and 
replacing two of them produces the 0(1040). The vector mesons empirically 
obey an equal spacing rule. This leads to a simple linear interpolation 
formula as a function of the light quark mass m: 

m v = muJ + (m* - muj ) ™\™ phys (29) 

171$ TYlphys 

The chiral limit is found with m = 0, resulting in 

m chiral = 783 - (1040 - 783) = 770 (30) 

TTl s Wlphys 

This is a 2% shift in the mass. 

One can address the quark mass dependence of the vector coupling to 
nucleons again through the SU(3) breaking pattern plus data. The SU(3) 
breaking of the NNV couplings themselves is not known well enough exper- 
imentally to be of direct use. However an indirect method can be employed. 
There is a phenomenologically successful model for these coupling, called 
vector meson dominance, that related the nucleon's NNV coupling g v to the 
coupling of the vector meson to a photon f v , defined via 

<V\J fM \0>=c v f v m 2 v e^ (31) 

where is the electromagnetic current and c v are known Clebsch-Gordon 
coefficients. The VMD relation is simply g v = f~ x , and can be moti- 
vated/derived from a dispersion relation argument for the electromagnetic 
current of the nucleon. Phenomenologically VMD works very well. Since 
there is experimental evidence on the SU(3) breaking of /„, we can use that 
plus VMD as an estimate of the quark mass dependence of g v . I have per- 
formed a fit to the data involving the neutral vector mesons p, u>, <fi again 
assuming a linear interpolation. The result is a dependence of the form 

9v = 9v0 [1 + b gv (ml - m 2 ph )] (32) 

with b gv = 0.57 ± 0.22 GeV -2 . The error bar comes from using different 
data in the evaluation of the SU(3)breaking. The central estimate is then 

9v \ chiral 



9v \ phy steal 



= 0.99 (33) 
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We can combine the vector mass and coupling to determine the shift in 
the parameter vector coupling Gy = gy/my. This yields 

Wt ch = Gv \cMral = q ggg ± Q.020 (34) 
*~*V \physical 

The variation of both the mass and the coupling constant were mild com- 
pared to the dependence that we will estimated for the scalar channel. More- 
over, both the mass in the denominator and the coupling constant in the 
numerator decreased, leading to a very small net effect in the ratio. The 
error bar quoted is a generous estimate. In any case, the uncertainty in this 
coupling is significantly smaller that that in the scalar coupling. 



8 Estimate of the changes in nuclear binding en- 
ergy 

Our results clearly have the strength of the attractive scalar coupling being 
significantly increased in the chiral limit. As I argued in Sec. II, in effective 
field theory the most consistent procedure would be to treat the long range 
components of the two pion exchange dynamically, most likely with a form 
of a cutoff, and then add the contact interactions to account for the rest 
of the interaction. This has not yet been done for heavy nuclei. However, 
the use of the contact interactions in point-coupling calculations of binding 
provide an estimate of the shift in the binding that is appropriate for the 
present state of the art. 

Not all possible contact interactions play a significant role in nuclear 
binding. The dominant ingredients in the binding of heavy nuclei have been 
elucidated by Furnstahl, Serot and co-workers 39 , 40 . The results can be 
extracted from Fig. 1 and Fig. 2 of jSH]- It is clear that the dominant 
effects are the scalar and vector contributions, to the leading power of the 
density. Other interactions play reduced roles, although for a complete 
understanding of the binding about a half-dozen contact interactions are 
required. I will consider only the dominant isoscalar-scalar and isoscalar- 
vector interactions, and in practice it is only the scalar coupling that has a 
significant dependence on the pion mass. 

Using Ref. [HS] , one can read off the effects of the different contact terms 7 
I parameterize the results in terms of the strengths of the contact interac- 

7 1 thank Dick Furnstahl and Brian Serot for assistance in understanding these numbers. 
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tions, normalized to their physical values, defining 

G s 



Vs 



Gs {physical 



W = (35) 

\physical 

The contributions to the binding energy numbers for 16 O (in MeV) are 

B.E. 



A 



-82t]s + Ur] V + 30 (36) 



The first two terms are the effects of the scalar and vector isoscalar inter- 
actions. The third term is the sum of 4 other smaller contributions to the 
binding energy and kinetic energy contributions. There is in addition the 
coulomb energy and a small center of mass correction. For 208 Pb, the result 
is 

~ -104775 + 57^ + 36 (37) 

The results of these calculations can be generalized to other nuclei by a 
parameterization that resembles the semi-empirical mass formula. For local 
interactions, because the nuclear density is nearly constant in the central 
region one expects that the binding energy will have a dependence volume, 
which in turn is proportional to the number of particles r 3 ~ A, and that 
interactions that occur near the nuclear surface would have a modified result 
proportional to the number of nucleons near the surface, r 2 ~ A 2 / 3 . This 
suggests that binding effects can be parameterized in terms of behavior in 
A and in A 2 / 3 . Using the results for nuclear matter and for specific nuclei, 
we find a good fit of the form 

B E 97 57 

— = -(120 - -^j^VS + (67 - -^ygOw + residual terms (38) 
Our results from the previous sections can be summarized as 

( Jr S\chiral -, - , n -. n 

Vs,ch = = 1.37 + 0.10 

[physical 

W:Ch = = 0.995 (39) 

\ physical 



These numbers produce 

B.E. 



A 



\chiral = -38 ± 8 MeV (40) 
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for O, compared to the physical result of —8 MeV. For Pb the corre- 
sponding results are — 49±10 MeV in the chiral limit compared to —11 MeV 
in the physical case. Finally for the general paramaterization, the results 
are suggest a shift in the binding energy 

A^Ural = "(44 - J^) (41) 

These shifts are larger than we would naively expect. There are two in- 
gredients in generating this magnitude. First is the large shift in the scalar 
contact interaction Gs- I have explained at length above why this large shift 
occurs and why it should be considered natural. The other ingredient is that 
there is a strong cancelation between the scalar and vector terms in the ex- 
pression for the binding energy. The usual 10 MeV/nucleon binding energy 
is the difference of two significantly larger numbers. This is understandable 
from the meson exchange potential description. Scalar exchange provides a 
strong attractive potential. Vector meson exchange provides the repulsive 
short range interaction and tends to oppose binding. The fact that the 10 
MeV of binding energy is far below all other hadronic scales of QCD comes 
from the competition of these two opposing effects. When we go to the chiral 
limit and make the attractive interaction significantly stronger while leaving 
the repulsive one unchanged, the near cancelation between these competing 
effects becomes less pronounced and the percentage shift is large. 



9 The chiral slope parameter 

Another way to present the results of this calculation is as the slope in the 
scalar coupling as one deviates from the chiral limit. The leading terms in 
an expansion around the chiral limit are given by 

Gs = G S0 (l + d s ml)+F s ml + ... (42) 

where Gso is the result in the chiral limit, and the non-analytic term is [20] 



f *=S}^-^> (43) 

The effect of the F$ term is relatively small, i.e. Fsm^/Gso ~ 0.09 . 

It should be pointed out that the dispersion relation does not exactly 
reproduce the correct m^. dependence of the chiral coupling. The cubic mass 
dependence arises in the dispersive calculation from the threshold behavior 
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of the chiral amplitude. Use of the chiral representation for the spectral 
function given in Eq. ^] leads to a cubic term that is 



I have checked both the spectral function and the real part of the loop 
calculation. This has nothing to do with the unitarization procedure and 
cannot be changed by higher order corrections to the threshold behavior. 
While the effect is numerically small, it represents a puzzle for which I do 
not understand the resolution 8 . It is fortunate that this nonanalytic term 
is numerically small and the effect of the disagreement has only a 2% effect 
on the determination of the slope. 

I have extracted the chiral slope by the evaluation of the dispersion 
integral at several values of the pion mass and performing a fit using terms 
in and m\ as well as the known nonanalytic term in m^. The result of 
the fit yields the slope 

4, = -(lT±6)feV-' = - a31± , a08 (45) 

The negative sign is indicative of the the fact that the scalar coupling is 
smaller in the physical case than it is in the chiral limit. This slope is 
relatively large for the reasons discussed above. Again, because the ratios of 
scalar couplings are predicted better than the absolute values, the parameter 
ds is much better determined in this calculation than D$ = dsGso- As seen 
in the previous section, the slope of the combination of the two isoscalar 
effects in the overall nuclear central force is yet larger because of the partial 
cancelation of the scalar and vector channels. 



10 Summary and discussion 

I have used the chiral results for the low energy behavior of the dispersion 
relation for the scalar contact interaction, supplemented with a representa- 
tion for the Omnes function, in order to calculate the strength of the central 
nuclear force in the chiral limit. In the results, the largest modifications 
are seen to come from the low energy amplitudes, which extend down to 
zero energy in the chiral case. While there are some uncertainties in the 
final result, most notably from the uncertainty in the chiral extrapolation of 

8 I thank E. Epelbaum, U. Meifiner, N. Kaiser and Barry Holstein for discussions about 
this issue 
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gAi the key ingredients in the calculation appear clear. We also understand 
qualitatively why the effect is relatively large. 

Since this has been an attempt to calculate the full scalar coupling, I 
have had to model the high energy contributions. While I have argued 
that this modeling has not introduced large effects in the final results, there 
is also a way to use this method with the full rigor of chiral perturbation 
theory. To do this, one would use the present calculation dynamically and in 
addition introduce a residual contact interaction to account for high energy 
effects which have been misrepresented in the present calculation. This extra 
parameter would presumably be small and would not be expected to have 
much mass variation. This variant of the Wilsonian scheme described in Sec. 
2 would then be a way to implement the effective field theory treatment while 
capturing the main dynamical results of the present calculation. 

There exists in the literature several other works which discuss the nu- 
clear interaction in the chiral limit [3J EJ within the context of chiral per- 
turbation theory. The present work is different because it is a dynamical 
attempt to calculate the mass variation in the scalar channel. Previous work 
has been both more general because this variation was parameterized a chi- 
ral coefficient, and also more limited because this coefficient could only be 
guessed at. In practice the mass variation found in this paper is larger than 
the estimates of jUEj, and has a well-determined sign. Another difference 
is that the previous works focussed on few-nucleon systems. This requires 
the addition of single pion exchange, which obviously also has a significant 
change in the chiral limit. My contribution here is only on the scalar sector 
- hence I have limited my comments on binding to the heavy nuclear case 
for which one pion exchange is not important. However, it is clear that 
the results of this paper would push the binding of few-nucleon systems in 
the direction of greater binding. It would be interesting to revisit the few- 
nucleon case with either the parameter dg calculated in the present work 
or with a full Wilsonian treatment dynamically including the effects of low 
energy two-pion exchange. 
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